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impinging  upon  spacecraft  surfaces.  Recent  experimental  evidence  from  the  SCATHA 
te  has  shown  that  charging  correlates  well  with  electrons  of  energies  greater  than 

In  this  paper  it  is  shown  that  the  SCATHA  observations  are  consistent  with  the 
if  charging  in  which  a  satellite  is  immersed  in  a  Maxwellian  plasma,  particle  collec- 
orbit  limited,  and  dominant  surface  effects  are  the  emission  of  secondary  and 
ttered  electrons.  The  energy  dependence  of  the  secondary  yield  for  multi-kilovolt  j 

t  electrons  determines  the  charging  threshold.  In  the  past,  inadequate  represen- 
of  the  secondary  yield  have  led  experimenters  to  question  the  validity  of  the 
g  model.  The  accuracy  of  the  secondary  electron  yield  formulation  based  on  electron 
g  power,  such  as  the  one  in  the  NASCAP  and  POLAR  charging  codes,  gives  good 
»nt  with  the  SCATHA  results.  A  Maxwellian  representation  of  the  magnetospheric  plasma 

c!?22sin9  effective. temperatures  and  densities  that  minimize  the  error  in 
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predictions  using  the  analytical  models  are  consistent. 
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1.  Introduction 


The  first  reports  of  large,  negative  spacecraft  surface  potentials  in  geosyn¬ 
chronous  orbit  were  of  eclipse  charging  (Deforest  1072;  Whipple  1981).  The 
charging  was  attributed  to  high  fluxes  of  keV  electrons  in  the  absence  of  photoem¬ 
ission,  and  explained  by  the  requirement  of  flux  balance  between  incoming  elec¬ 
trons,  backscattered  electrons,  secondary  electrons,  photoelectrons,  and  incoming 
ions.  Early  efforts  to  correlate  spacecraft  potentials  with  parameters  describing 
the  incident  electron  population  used  either  constants  or  fits  to  low  energy  (0-1000 
eV)  data  to  calculate  the  backscatter  and  secondary  yield  coefficients  of  the  space¬ 
craft  surface  materials.  (Garrett  and  Rubin,  1978;  Garrett  and  DcForest,  1979; 
Garrett  et  al.,  1980).  These  efforts  could  not  consistently  predict  charging  levels. 

Subsequently,  charging  events  in  excess  of  one  kilovolt  negative  were  meas¬ 
ured  in  sunlight.  Sunlight  charging  is  more  complex  than  eclipse  charging 
because  of  the  lack  of  spatial  symmetry.  It  is  clear  that  additional  processes  are 
at  work,  since  first  order  calculations  of  the  ambient,  secondary,  and  photo- 
currents  produce  net  postive  currents  to  a  magnetospheric  satellite.  Direct  evi¬ 
dence  from  measured  particle  spectra,  and  inferred  evidence  from  modeling  efforts, 
suggest  that  the  formation  of  potential  barriers  around  the  satellite  suppress  pho¬ 
toemission  when  the  electron  environment  is  sufficiently  intense  ( Whipple , 
1976a, b;  Oleen  et  al.,  1981;  Mandell  et  al.,  1978).  Generally,  a  given  environment 
is  expected  to  charge  a  satellite  about  an  order  of  magnitude  more  in  eclipse  than 
in  sunlight.  ( Garrett  (1980)  states  that  "Potentials  as  high  as  -2000  V  in  sunlight 
and  -20,000  V  in  eclipse  have  been  observed  on  ATS-6."  See  also  Purvis  (1983).) 
This  is  because,  although  current  balance  to  a  shadowed  insulator  is  the  same  as 
that  to  an  eclipsed  spacecraft,  the  charging  rate  is  typically  three  to  four  orders  of 
magnitude  slower  since  the  insulator  must  differentially  charge  relative  to  under¬ 
lying  grounded  metal.  It  follows  that  the  maximum  potential  achieved  in  a 


1 


daylight  charging  errant  it  governed  by  the  duration  of  the  charging  environment, 
and  ia  usually  far  below  the  "equilibrium"  (steady  state)  value  which  might  be 
predicted.  Olsen  and  Purvis  (IMS)  discuss  the  charging  dynamics  for  daylight 
and  eclipse  charging  events  observed  on  the  ATS-5  and  ATS-6  spacecraft. 

A  number  of  studies  have  attempted  to  correlate  charging  with  specific 
features  of  observed  electron  distributions.  Reagan  it  of.  (1M1)  found,  for  the 
SCATHA  satellite,  that  "surface  potentials  ...  were  determined  primarily  by  ... 
electrons  in  the  energy  range  <  SO  keV  and  ...  that  surface  charging  occurred 
when  the  spectrum  hardened."  Olsen  (IMS)  found,  in  studies  of  ATS-6  and 
SCATHA  data,  that  when  the  electron  spectrum  hardened  his  "count  rate"  com¬ 
monly  exhibited  a  sharp  drop  just  above  its  maximum,  and  that  charging 
occurred  when  this  drop  was  at  an  energy  exceeding  about  15  keV. 

Recently,  a  comprehensive  study  of  high-level  sunlight  charging  events  on 
the  SCATHA  satellite  by  Mullen  et  at.  (1080)  showed  that  the  measured 
spacecraft- to-plasma  potential  difference  is  directly  proportional  to  the  intensity 
of  the  ambient  electron  flux  greater  than  30  keV,  and,  additionally,  that  the 
spacecraft  potential  is  not  consistently  related  to  the  electron  flux  below  this 
energy,  despite  the  fact  that  the  low  energy  flux  is  commonly  an  order  of  magni¬ 
tude  higher  than  the  high  energy  flux.  (Note  that  the  spectra  exhibited  by  Olsen 
(1983)  also  show  significant  counts  at  30  keV  and  above,  despite  a  sharp  drop  in 
the  10-20  keV  range.)  A  similar  relationship  between  high  energy  electron  flux 
and  spacecraft  potential  in  excess  of  -100  volts  was  found  for  low  altitude  eclipse 
charging  of  DMSP  satellites  in  polar  orbit  (Gussenhoven  et  al 1985).  These 
results  were  explained  in  terms  of  the  concept  of  a  critical  energy  below  which  the 
electrons  do  not  contribute  to  charging  because  they  are  self-balanced  by  their 
own  secondary  and  backscattered  production.  The  concept  of  critical  energy  was 
quantitatively  developed  for  incident  Maxwellian  electron  distributions  by  Lai  et 
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al.  (1083). 

The  net  current  density  to  *  satellite,  associated  with  incident  elec¬ 
trons,  is  found  by  subtraetinf  backscatter  and  secondary  electrons  from  the 
incident  electron  distribution  and  integrating  over  energy.  If  J™*  is  negative,  the 
spacecraft  will  achieve  a  negative  potential  such  that  is  balanced  by  incident 
ions  and  ion-produced  secondary  electrons. 

Previous  studies  by  several  authors  have  shown  that,  for  a  Maxwellian 
plasma,  charging  results  only  when  the  electron  environment  exceeds  some  thres¬ 
hold  temperature  ( Stannard  et  a/.,  1981;  Laframboise  and  Kamitaumc,  1083). 
This  is  because  secondary  yields  typically  are  greatest  for  electrons  with  energies 
below  about  1  keV  and  are  less  than  unity  for  electrons  with  energies  above  a  few 
keV.  Below  the  threshold,  the  integral  secondary  and  backscatter  yields  exceed 
the  incident  electron  flux,  so  net  charging  cannot  occur. 

Both  the  study  by  Laframboiat  and  Kamitauma  (1083),  and  that  by  Lai  et  al. 
(198S),  indicated  charging  threshold  temperatures  of  a  few  keV,  and  that  charging 
was  caused  by  electrons  with  energies  on  the  order  of  10  keV.  These  energies  are 
about  a  factor  of  three  lower  than  the  threshold  energy  inferred  from  the  SCA- 
T3A  data  (Afutien  et  al .,  1988).  As  shown  below,  this  discrepancy  occurs  because 
the  formulation  of  secondary  yields  used  by  these  authors  was  inaccurate  and 
predicted  unphysically  small  secondary  yields  for  electrons  in  the  relevant  energy 
range  (5-50  keV). 

Lai  et  al.  (1983)  and  Laframboiat  and  Kamitauma  (1983)  used  an  analytical 
expression  for  the  secondary  production  given  by  Sandera  and  Inouye  (1979) 
derived  from  the  work  of  Stemglaaa  (1954).  The  secondary  production  predicted 
is  accurate  primarily  for  electrons  of  low  energy  (i.e.,  up  to  a  few  keV)  but  is 
essentially  zero  for  electron  energies  greater  than  10  keV.  Theory  ( Alig  and 
Bloom ,  1978)  and  experiment  ( Kanter ,  1961)  confirm  that  secondary  yield  is 


proportional  to  stopping  power  for  high  energy  incident  electrons.  Using  tec  on* 
dery  formulations  which  ere  more  correct  at  high  energy,  such  as  those  from  NAS- 
CAP  (Kcts  ft  el.,  1977;  MandfU  ft  el.,  1984),  or  that  suggested  by  Burke  (Burke, 
1980),  better  agreement  with  the  SCATHA  obeerratione  is  obtained. 

la  the  following  we  formulate  a  more  accurate  expression  for  seconder/  pro¬ 
duction  from  the  high  energy  electron  population  measured  at  geosynchronous 
orbit.  We  propose  a  method  for  selecting  Maxwellian  parameters  from  arbitrary 
distribution  functions  which  best  calculates  electron  charging  current  densities. 
The  observed  strong  correlation  of  <»Wgi»g  with  the  flux  of  electrons  with  energy 
greater  than  30  keV  and  the  non-existent  (if  not  negative)  correlation  with  elec¬ 
trons  below  30  keV  is  easily  demonstrated  when  the  electron  environment  is 
modeled  with  distributions  having  constant  density  and  various  temperatures. 

3.  Theory 

The  theory  developed  in  this  section  is  for  a  spherical  object  in  an  isotropic 
Maxwellian  plasma  under  eclipsed  conditions.  However,  the  conclusions  drawn 
are  reasonably  general.  AU  quantities  are  in  SI  units,  with  particle  energies  and 
plasma  temperatures  in  eV,  unless  otherwise  stated. 

The  time-dependent  charging  of  an  object  immersed  in  a  plasma  whose  total 
current  density  to  the  object  is  J(Qt  can  be  described  by  the  equation 

dV/dt  -  R  Jtot  /  «0  , 

where  V  is  the  object  to  plasma  potential  difference,  R  is  the  effective  radius  of 
the  object,  and  is  the  permittivity  of  free  space.  When  in  eclipse,  the  current 
density  consists  of 

Jtot  “  ^el  "**  -Hon  ***  ^sec  +  ^back  +  ^prosec  * 

Here  j^  and  j.QQ  are  the  incident  electron  and  ion  current  densities  from  the 
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ambient  environment,  j##c  and  ere  the  electron  secondary  and  baekacattered 
current  densities  produced  by  the  incident  electron  current  density,  and  j  is 

pros«v 

the  electron  secondary  current  density  produced  by  the  incident  ion  current  den¬ 
sity.  When  the  Incident  election  population  can  be  described  by  a  Maxwellian 
distribution  with  temperature  #,  the  secondary  Q  )  and  bachscatter  (j^^) 
current  densities  can  be  written  as  integrals  over  the  incident  spectrum: 

WJ«1  "  ‘'Vo"  E  «P(-E/»)  V(S)  d£ 

"  E  «5>(‘e/*)  S(E)  dE  . 

In  these  equations  Y(E)  and  B(E)  are  yield  functions  for  secondary  and  backseat- 
tered  production,  respectively,  taking  into  account  the  angular  distribution  of  the 
incident  electrons. 

2.1.  Secondary  Emission  Yields 

Typical  yield  of  secondary  electrons  as  a  function  of  energy  is  shown  in 
figure  1.  The  yield  is  proportional  to  the  energy  deposited  by  the  incident  elec¬ 
tron  in  the  top  20-100  angstroms  of  the  material.  For  low  energy  electrons  (wi*h 
range  comparable  to  or  less  than  this  distance)  the  yield  is  proportional  to  the 
particle  energy.  At  high  energy  the  yield  u  proportional  to  the  stopping  power. 
This  behavior  has  both  a  strong  theoretical  basis  and  has  been  verified  experimen¬ 
tally  ( Kanter ,  1961).  The  maximum  yield  occurs  in  the  transition  between  these 
two  regimes,  which  is  at  a  few  hundred  volts  for  most  materials.  For  spacecraft 
charging  purposes  it  is  important  to  have  a  reasonable  estimate  of  secondary 
yields  for  5-50  keV  primary  electrons.  For  most  materials  the  electron  range  is 
well-represented  in  this  energy  regime  by  a  power  law: 

Range  —  Constant  x  Energyp 


with  1.5<p<2.0.  Since 


Stopping  Power  —  [d(Range)/d(Euergy))‘l 

« 

wt  npwt  secondary  yield  to  M  off  inversely  (or  a  lower)  with  reaped  to  the 
energy  of  the  primary  electron: 

Y(E)  -  l*/SKtt]1‘p 

where  E«t  “  the  energy  at  which  the  aecondaxy  yield  extrapolatee  to  unity.  ] 

(This  is  equivalent  to  the  form  suggested  by  Burke  (1980),  who  used  p«1.725. 

Taking  into  account  that  Burke's  data  is  for  normal  incidence  whereas  we  are 
interested  in  isotropically  incident  electrons,  the  relation  between  Burke's  con¬ 
stant  K  and  B^  is 

B«rt  “  («)l  JTM  |k.V|. 

The  Appendix  suggests  values  of  p  and  B  ,  for  several  materials.)  Proposed 

ext  , 

simplified  formulas  (Sanders  and  /nouye,  1979;  Sternglass,  1954)  having  exponen-  4 

tial  falloffs,  while  easy  to  use  analytically,  will  invariably  underestimate  the  yield  I 

of  secondary  electrons  for  environments  capable  of  charging  spacecraft.  For  most  ^ 

spacecraft  materials  (which  have  low  atomic  number)  backscatter  yields  increase  j 

monotonically  with  energy  from  sero  to  a  few  tenths.  J 

| 

Following  Lai  et  el.  (1983),  E  is  defined  as  the  energy  below  which  ' 

upper 

electrons  do  no  charging,  i.e., 

t 

/0®«pp«  E  exp(-E/0)  (l-Y(E)-B(E))  dE  -  0  j 

1 

The  parameter  &UDOtt.  is  a  function  of  9 ,  being  undefined  for  9  below  the  charging  ; 

threshold,  infinite  at  threshold,  and  reaching  a  finite  limit  as  9  ->  00.  Assuming  "  \ 

\ 

the  high-energy  formula  given  above  for  Y(E)  and  a  constant  backscatter  J 

coefficient,  B,  this  limiting  value  is  < 

-  E«*  i!*-P)(l-B)/2]1/(1-P>  : 


For  Gold  subject  to  isotropically  incident  electron*,  we  suggest  (based  on  values 
developed  for  NAS  CAP,  Katz  at  of.,  1077;  MandaU  at  o 1084), 
the  values  E^  •*  4.6  keV,  B  -  .68,  end  p  »  1.73.  We  then  find  EUpper  “  36 
keV.  Using  the  fall  secondary  emission  formulation  found  in  NASCAP  gives  a 
•lightly  lower  value  of  31.2  keV.  This  contrasts  with  the  11.1  keV  value  reported 
by  Lai  et  at .  (1083).  The  threshold  Maxwellian  temperature  for  charging  U 
found  by  NASCAP  to  be  14.0  keV,  compared  with  4.0  keV  from  Lai  et  al. 
(1083).  The  low  values  reported  by  Lai  are  due  principally  to  the  inadequacy  of 
the  secondary  formulation  used.  The  backscatter  coefficient  used  here  is  also 
somewhat  higher  than  used  by  Lai  et  al.  (1083). 

For  the  purpose  of  comparison  with  data  from  the  SCATHA  satellite,  the 
electron  emission  properties  are  taken  to  be  those  of  soler  cell  rover  slip  model 
developed  for  NASCAP  ( MandaU  at  el,  (1084))  since  most  of  the  exterior  surface 
of  the  spacecraft  is  covered  with  solar  cells.  (The  results  would  change  little  for 
most  typical  spacecraft  coverings  such  as  teflon  beta  cloth.)  Here  the  properties 
E^  »  4.6  keV,  B— .32,  and  p  —  1.73,  used  in  the  formula  above,  give  a  value 
Eupp«r(°°)  "  14*3  keV.  The  dashed  curve  in  figure  2  shows  the  dependence  of 
Eupper  (c*lcul*t*4  using  the  full  NASCAP  treatment  for  secondary  and  back¬ 
scatter  yields)  with  the  temperature  of  the  ambient  environment.  The  limiting 

value  of  E  by  this  method  is  15.1  keV,  and  the  threshold  temperature  is  6.8 
upper 

keV.  These  are  much  lower  than  the  corresponding  values  for  Gold,  which  is  & 
high  atomic  number  material.  The  average  energy  of  the  charging  electrons,  that 
is  those  above  E  weighted  by  their  charging  effectiveness  is  given  by 


J  09  E2  exp(-E/f)  [1-Y(E)-B(E)1  dE 

E  •»  — — — — — — — — — — — — — 

charging 

/  "  E  exp(-E/#)  (l-Y(E)-B(E)]  dE 

The  solid  curve  in  figure  2  shows  the  dependence  of  this  value  on  temperature. 
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Unlike  S _ ,  the  iwin  value  has  an  minimum  near  the  threshold  tempera* 

upper 

ture.  For  tl±e  SCATHA  materials  the  minimum  value  is  35  keV  at  a  temperature 

of  0  keV.  The  flux  of  electrons  with  energies  greater  than  E _ and  within  9  of 

upper 

the  Elqhargiag  contribute  most  of  the  net  electron  current  to  the  surface  and 
should  be  a  reliable  indicator  of  charging.  The  choice  of  is  not  unique; 

the  median  value  above  E  would  suffice.  The  basic  concept  is  just  that 
charging  is  driven  by  electrons  above  and  within  a  few  9  of  E  . 

The  quantities  Eupp<r  and  ErW^ng  relate  to  the  findings  of  Olsen  (1983) 
and  Mullen  et  at.  (1986)  in  the  following  manner.  A  given  environment  will  not 
charge  a  material  if  the  actual  electron  spectrum  drops  significantly  below  its 
Maxwellian  approximation  at  an  energy  much  less  than  the  material's  E 

upper' 

current  balance  would  then  be  dominated  by  the  secondary  electrons  produced  by 

the  low  energy  part  of  the  spectrum.  Since  E„_ _ was  calculated  above  to  be 

upper 

15.1  keV,  this  explains  the  finding  of  Olsen  (1983)  that  the  requirement  for  charg¬ 
ing  is  that  such  a  drop  occur  at  15  keV  or  higher.  Given  this,  one  then  expects 
the  spacecraft  potential  to  correlate  with  the  net  charging  current.  By  definition, 
the  net  charging  current  consists  of  that  portion  of  the  spectrum  with  energy 
above  Eupper,  and  is  naturally  correlated  with  the  differential  flux  at  E^^g^g- 
This  explains  tbe  finding  of  Mullen  et  al.  (1986)  that  spacecraft  potential  is  corre¬ 
lated  with  the  flux  at  30  keV. 

3.3.  Single  Maxwellian  Representation  of  Environment 

The  prece'/jig  discussion  of  net  current  density  has  been  for  an  assumed 
Maxwellian  distribution  of  electrons.  Data  indicates  that  charging  environments 
are  not  well  fit  by  a  single  Maxwellian;  frequently  two  or  more  Maxwellian 
descriptions  are  far  superior.  The  net  electron  charging  current  density  is  in  gen¬ 
eral  the  integral  of  the  first  moment  of  the  velocity  distribution  function  with  a 


weighting  function  equal  to  one  minus  the  secondary  plus  backseat  ter  yields. 
That  ie, 

"  /.- E  w  W(E)  dE» 

where  W(E)  —  (l-Y(E)-B(E)].  For  a  Maxwellian  distribution  the  function  f(E)  is 

f(E)  -  N  e  f2  (etf/2xmj1/2  exp(-E/0) 

where  N  is  the  number  density,  m  and  e  are  the  mass  and  the  charge  of  the  elec- 
tron,  respectively.  The  ideal  Maxwellian  fit  (for  charging  purposes)  would  have 
the  property  of  exactly  matching  However,  since  calculating  Jg  requires  a 

computer  program  even  for  a  Maxwellian,  a  simpler  recipe  is  highly  desirable. 


In  the  energy  range  of  interest  (a  few  to  fifty  keV)  the  backscatter  yield, 
B(E),  depends  weakly  on  energy,  while  the  secondary  electron  yield  Y(E)  is  rea¬ 
sonably  fit  by  a  power  law  E1"15.  Thus  we  choose  a  density  and  temperature 
which  match  the  incident  electron  flux  and  the  flux-weighted  mean  of  E1-p. 
Matching  the  incident  flux  assures  a  reasonably  close  value  for  the  weakly  energy 
dependent  (l-B(E)j  portion  of  the  integral  for  J  ,  while  the  other  condition  gives  a 
good  match  for  the  mean  secondary  electron  yield.  For  a  Double  Maxwellian  with 
densities  N^,  Ng,  and  temperatures  9 0g,  respectively,  we  define 

FA-NA*Al/2^B-NB»Bl/J- 


Then,  it  can  be  shown  that  the  single  Maxwellian  which  fits  our  criteria  has  tem¬ 
perature  and  density  given  by 


»l'P  -  (FA/(FA+FB))«A1-p  +  (FB/(FA+FB))*Bl-p 
N  -  [FA+FB]f1/2 . 


We  have  applied  this  scheme  to  fifty-three  Double  Maxwellian  environments 
measured  by  the  SCATHA  spacecraft,  as  quoted  by  Sehnuelle  et  al.  (1981).  We 
used  the  value  p«1.73,  appropriate  to  the  NASCAP  model  for  solar  cell 


cowslips,  as  wall  ma  many  other  common  spacecraft  materials.  The  full  NAS- 
CAP  formulations  for  secondary  and  bacicscatter  yield  were  then  used  to  calculate 
J^et  for  each  Double  Maxwellian  environment  and  its  single  Maxwellian  represen¬ 
tation.  Defining  the  error  as  the  difference  in  divided  by  the  incident  elec¬ 
tron  flux;  we  found  the  root-mean-square  error  to  be  3.7  percent,  with  the  max¬ 
imum  error  9.2  percent.  The  larger  errors  were  for  environments  below  the  thres¬ 
hold  for  charging.  The  five  environments  capable  of  charging  this  material  had  a 
root-mean-square  error  of  0.6  percent. 

By  this  procedure  the  equivalent  Maxwellian  temperature  for  charging  is 
material  dependent  through  the  range  exponent  p.  However,  since  the  variation 
in  p  is  not  great,  this  material  dependence  is  fairly  slight.  Also,  since  secondary 
electrons  produced  by  ions  are  proportional  to  the  ion  velocity,  the  value  p»l/2 
should  be  used  to  compute  the  equivalent  Maxwellian  temperature  for  an  ion 
spectrum. 

2.3.  Application  to  Spacecraft  Charging 

We  will  now  show  that  the  above  considerations  predict  a  positive  correlation 
between  spacecraft  potential  and  measured  electron  flux  with  energy  in  excess  cf 
30  keV,  and  a  null  (or  even  negative)  correlation  with  electron  flux  below  30  keV 
(or,  equivalently,  total  electron  flux). 

Consider  tbe  charging  response  of  an  idealized  satellite  in  a  neutral  isother¬ 
mal  Maxwellian  hydrogen  plasma.  The  satellite  is  assumed  to  be  a  spherical 
probe  that  collects  ion  and  electron  currents  according  to  orbit  limited  theory. 
That  ia,  net  current  densities  resulting  from  the  incident  electron  and  ion  distribu¬ 
tions  e  .*e 


Tnet 

Jel 


E  f  (E-V)  W(E)  dE, 
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I 

VMLr#7 


Jj  -  ;.v*  E  fj(E+V)  Wj(E)  dE 

where  V  is  the  spacecraft  potential  (assumed  negative),  Wj(E)  is  unity  plus  the 
secondary  electron  emission  coefficient  for  incident  ions,  and 

f0(s)  -  •  r2  «*p(-E/») 

The  equilibrium  potential  is  found  by  integrating  in  time  the  bask  charging  equa¬ 
tion  until  a  steady  value  is  obtained.  (The  time  scale  for  approach  to  steady  state 
is  milliseconds  for  typical  magnetospheric  plasmas.)  The  equilibrium  floating 
potential  is  a  monotonic  function  of  the  plasma  temperature. 

-3 

For  a  fixed  plasma  density  of  1  cm  the  equilibrium  satellite  potential  was 
calculated  for  plasma  temperatures  ranging  from  8  to  20  keV  in  steps  of  one  kilo¬ 
volt.  Also  calculated  was  the  incident  electron  flux  in  the  two  energy  ranges  of  the 
instrument  on  SCATHA  reported  in  Mullen  et  al.  (1986).  The  low  energy  chan¬ 
nel  was  taken  to  be  all  electrons  between  50  eV  and  30  keV  the  high  channel  was 
from  30  keV  to  400  keV. 

In  Figure  3  the  circles  indicate  spacecraft  potential  as  a  function  of  the  high 
energy  flux,  the  plus  signs  the  potential  vs  the  low  energy  flux.  Each  pair  of  cir¬ 
cles  and  crosses  at  a  single  potential  are  results  irom  a  single  temperature.  The 
obvious  positive  correlation  of  the  potential  with  high  energy  flux  and  the  small, 
negative  correlation  with  the  low  energy  flux  are  in  good  agreement  with  the 
SCATHA  data  from  Mullen  et  al.  (1986).  (The  scale  of  the  potential  here  is  an 
order  of  magnitude  greater,  because  tbia  analysis  was  performed  for  eclipse  charg¬ 
ing,  while  the  published  observations  are  for  sunlight  charging.)  The  incident 
electron  current  in  all  cases  of  charging  is  balanced  by  a  combination  of  incident 
ions  and  secondary  and  backscattered  electrons. 
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3.  Diinutioa 

The  observation  that  charging  correlates  directly  with  the  flux  of  electrons 
with  energies  greater  than  30  keV  is  completely  consistent  with  the  description  of 
charging  as  the  balancing  of  incident  electron  currants  by  secondary  emission, 
baekscatter  and  ion  currents.  Previous  analysis  by  Lai  et  at.  (1083)  had  assumed 
the  material  properties  were  that  of  gold,  a  highly  emissive  material,  but  still 
came  up  with  lower  than  the  observed  energies.  The  present  work  is  in  much 
better  agreement  with  experiment  for  two  reasons.  First,  and  certainly  most 
important,  the  formulation  of  the  secondary  yield  as  a  function  of  energy  used 
here  is  much  more  accurate  for  incident  electron  energies  greater  than  a  few  keV. 
Second,  the  definition  of  the  mean  energy  of  the  charging  electrons  is  more 
relevant  than  that  of  the  lowest  energy  electron  to  contribute  to  charging,  since 
the  bulk  of  the  charging  current  is  near  the  mean.  The  fundamental  mechanisms 
leading  to  spacecraft  charging  are  the  same  for  both  works,  namely,  the  achieve¬ 
ment  of  current  balance.  However,  it  points  out  just  how  critical  is  the  accurate 
knowledge  of  both  material  properties  and  the  ambient  environment  for  both  high 
and  low  energies  if  modeling  efforts  are  to  give  good  spacecraft  potential  predic¬ 
tions.  While  this  paper  used  a  single  Maxwellian  representation  of  the  charging 
environment,  the  results  are  applicable  to  other  energy  distributions  since  charg¬ 
ing  involves  integrals  over  the  distribution  and  is  not  particularly  sensitive  to 
details  of  the  shape  of  the  distribution  function.  Tbs  prediction  of  the  energy 
channel  which  will  indicate  charging  will  vary  little  for  well  behaved  monotonic 
spectra.  The  good  agreement  with  SCATHA  results  from  having  a  good  represen¬ 
tation  of  the  material  properties  and  the  high  energy  electron  environment.  Use 
of  a  single  Maxwellian  to  represent  the  highly  variable  magnetospheric  environ¬ 
ment  is  convenient,  but  not  essential. 


Appendix 

Suggested  Constant*  for  Several  Materials 

The  table  below  gives,  for  several  materials,  suggested  values  for  the  range 
exponent  p,  the  energy,  E^,  at  which  secondary  yield  extrapolates  to  unity, 
approximate  backscatter  coefficient,  B  (evaluated  at  10  keV),  and  resultant  limit¬ 
ing  value  of  E _ from  the  formula 

upper 

Eupp«(oo)  “  Eext  ((3-p)(l-B)/2|1''<1-p) 

Values  from  Burke  (1980)  have  been  doubled  to  account  for  isotropic  incidence. 


Material 

P 

e^PmV! 

B 

E  [keVl 

upper1  J 

Aluminum 

1.78  (B) 

1.8  (B) 

.38(B) 

6 

Carbon  (Aquadag) 

1.55  (B) 

1.2  (B) 

.27(B) 

4 

Gold 

1.73  (B) 

4.8  (B) 

.64(B) 

35 

Kapton 

1.725(A) 

1.53(A) 

.25(B) 

4 

Lucite 

1.725(A) 

3.02(A) 

■14(A) 

7 

Magnesium 

1.75  (B) 

.7  (B) 

.35(B) 

2 

Mylar 

1.725(A) 

2.07(A) 

•14(A) 

5 

Nylon 

1.725(A) 

3.15(A) 

•14(A) 

7 

Polyethylene 

1.725(A) 

4.02(A) 

■14(A) 

9 

Polystyrene 

1.725(A) 

2.27(A) 

•14(A) 

5 

Polyvinylalcohoi 

1.725(A) 

3.51(A) 

•14(A) 

8 

Silver 

1.74  (B) 

4.0  (B) 

.55(B) 

22 

Si02 

1.86  (B) 

4.8  (B) 

.33(B) 

15 

SOLA 

1.73  (B) 

4.63(B) 

.33(B) 

15 

Teflon 

1.725(A) 

4.75(A) 

.29(B) 

14 

(A)  Burke  (1980) 

(B)  NASCAP  (Katz  et  al.,  1977;  Mandcll  et  al.,  1984) 
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Figure  1.  Secondary  yield  vs.  primary  electron  energy  for  solar  cell  coverslip 
material,  assuming  isotropic  incidence. 
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Fifure  2.  (dashed  curre)  and  (solid  curve)  as  functions  of 

assumed  Maxwellian  temperature 
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